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Abstract 

An  improved  zero-effort-miss  (ZEM)  and  zero-effort-velocity  (ZEV)  optimal  feedback  guidance  algorithm  is  proposed  in  this  paper 
for  Mars  powered  descent  phase,  which  has  taken  into  account  the  no-subsurface  constraints  when  designing  the  guidance  law.  Detailed 
descriptions  have  been  given  to  this  algorithm.  The  main  advantage  of  this  optimal  guidance  law  is  that  it  can  successfully  avoid  collision 
to  the  Martian  surface,  at  the  same  time  preserving  high  landing  precision  and  fuel  optimality.  For  future  Mars  landing  mission,  the 
thrust-limited  engine  is  most  likely  to  be  used.  We  propose  a  simulation  method  to  study  and  determine  the  optimal  final  time  which 
can  minimize  the  fuel  usage  while  ensuring  precision  landing  and  no-subsurface  constraints.  Monte  Carlos  are  run  to  test  the  robustness 
and  reliability  of  the  improved  and  original  guidance  algorithms  which  show  the  superiority  of  the  improved  guidance  law  in  preventing 
collision  to  the  surface  over  the  original  one. 

©2014  COSPAR.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Landing  on  Mars  has  been  a  hot  research  area  in  recent 
years.  Many  efforts  have  been  made  to  it,  such  as  the  Vik¬ 
ing,  MPF,  Spirit,  Opportunity,  Phoenix  and  MSL  mis¬ 
sions.  The  Mars  entry,  descent,  landing  (EDL)  mission 
roughly  includes  four  phases:  the  hypersonic  entry  phase, 
the  subsonic  parachute  descent  phase,  the  powered  descent 
phase  and  the  touchdown  phase  (Braun  and  Manning, 
2007).  The  purpose  of  EDL  mission  is  to  deliver  a  vehicle 
to  a  desired  place  on  Martian  surface,  which  requires 
pinpoint  landing.  The  Mars  entry  guidance  is  a  tough 
problem  due  to  the  thinness  and  high-uncertainties  of  Mar¬ 
tian  atmosphere.  The  purpose  of  the  entry  guidance  is  to 
drive  a  vehicle  to  a  desired  parachute  deployment  site. 
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The  parachute  is  deployed  and  this  parachute  descent 
phase  is  totally  uncontrolled.  Then  the  parachute  is  dis¬ 
carded  and  the  vehicle  enters  the  powered  decent  phase, 
which  will  be  discussed  in  this  paper  in  detail.  The  landing 
guidance  should  be  designed  to  satisfy  the  requirement  of 
precision  landing  and  fuel  optimality.  Many  guidance  laws 
and  optimization  methods  have  been  proposed  to  address 
this  problem.  There  are  now  many  optimization  softwares 
on  the  market,  such  as  GPOPS,  SNOPT,  DIDO,  TOM- 
LAB  and  others.  These  softwares,  which  are  based  on  var¬ 
ious  optimization  methods,  can  automatically  design  an 
optimal  trajectory  which  can  satisfy  the  motion  equations 
and  various  constraints  at  the  same  time. 

Convex  optimization  method  can  also  be  used  to  design 
an  optimal  trajectory  for  Mars  landing  problem  with  many 
constraints  (Acikmese  and  Ploen,  2007)  and  Blackmore 
et  al.  (2010)  improved  the  convex  optimization  method  to 
generate  the  minimum  landing  error  trajectory  when  no 
feasible  trajectory  to  the  target  exists.  However,  the 
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powered  decent  guidance  based  on  convex  optimization  is 
an  open  loop  guidance  which  cannot  be  robust  against  dis¬ 
turbance,  so  an  optimal  feedback  guidance  law  is  needed 
for  the  landing  mission  which  cannot  only  be  fuel-optimal, 
but  also  achieve  the  precision  landing  requirement. 

Ebrahimi  et  al.  (2008)  proposed  a  new  optimal  sliding 
mode  guidance  control  law  based  on  the  zero-effort-miss 
(ZEM)  and  zero-effort- velocity  (ZEV).  The  zero-effort- 
velocity(ZEV)  error  is  a  new  concept  proposed  by  Ebrah¬ 
imi  et  al.,  which  is  analogous  to  the  zero-effort-miss 
(ZEM)  distance.  ZEV  is  defined  as  the  velocity  error  at 
the  final  time  of  the  mission  if  no  commanded  acceleration 
is  added  to  the  vehicle  since  the  current  time.  Ebrahimi 
et  al.  integrated  this  optimal  guidance  law  with  sliding- 
mode  control  for  fixed-interval  propulsive  maneuvers  with 
terminal  velocity  constraint  (Ebrahimi  et  al.,  2008). 
Furfaro  et  al.  (2011)  later  modified  this  ZEM/ZEV  optimal 
guidance  law  for  a  lunar  precision  landing  mission  which 
has  shown  robust  performance  against  dispersions  by 
simulation.  But  besides  the  precision  landing  and  fuel- 
optimality  requirement,  the  landing  mission  also  requires 
no-subsurface  flight.  The  ZEM/ZEV  optimal  guidance 
law  does  not  account  for  this  constraint,  which  may 
cause  the  vehicle  to  crash  the  ground  following  this  guid¬ 
ance  law.  In  order  to  rule  out  this  possibility,  Guo  et  al. 
(2013)  proposed  a  waypoint-optimized  zero-effort-miss/ 
zero-effort-velocity  feedback  guidance  for  Mars  landing, 
which  can  autonomously  design  an  optimal  waypoint  using 
optimization  approach  during  the  powered  descent  phase  if 
a  subsurface  flight  is  predicted.  The  designed  waypoint  is 
regarded  as  an  intermediate  target  and  the  ZEM/ZEV 
algorithm  is  employed  to  reach  it,  and  then  reach  the 
desired  landing  site. 

The  subsurface  flight  is  most  likely  to  happen  if  the  land¬ 
ing  site  is  behind  the  vehicle’s  initial  position  of  powered 
descent  phase  which  is  a  rare  condition.  By  carefully 
designing  the  reference  trajectory,  this  possibility  can  be 
reduced  or  even  eliminated.  So  in  most  of  the  cases,  the 
landing  site  is  ahead  of  the  vehicle’s  initial  position. 
Besides,  when  the  subsurface  flight  is  predicted,  whether 
the  optimized  waypoint  can  be  found  in  a  short  period  of 
time  is  also  a  crucial  question. 

So  in  this  paper,  we  improve  the  zero-effort-miss  (ZEM) 
and  zero-effort-velocity  (ZEV)  optimal  guidance  law  by  tak¬ 
ing  the  no-subsurface  constraints  into  account  when  design¬ 
ing  the  guidance  law.  This  is  a  trade-off  between  the 
precision  landing,  fuel  optimality  and  no-subsurface 
requirements.  We  sacrifice  a  little  landing  precision  to  avoid 
the  possibility  of  collision  to  the  surface.  According  to  the 
improved  ZEM/ZEV  optimal  guidance  law,  the  thrusting 
force  can  be  automatically  changed  to  prevent  collision  to 
the  surface.  Detailed  descriptions  have  been  given  this  algo¬ 
rithm.  Besides,  a  new  experimental  study  has  been  proposed 
to  study  the  effect  of  different  final  time  on  the  performance 
of  the  optimal  feedback  guidance  regarding  the  fuel-opti¬ 
mal,  pinpoint  landing  and  no-subsurface  requirement  and 
the  optimal  final  time  is  determined  by  trial.  In  simulation 


study,  the  target  terminal  altitude  of  the  powered  descent 
phase  is  set  to  100  m,  which  is  rational  considering  the  alti¬ 
tude  needed  for  the  final  touchdown  phase.  And  the  pow¬ 
ered  descent  phase  of  the  MSL  mission  also  stops  at  100  m 
to  let  the  sky  crane  touchdown.  Monte  Carlos  are  run  to  test 
the  robustness  and  reliability  of  the  improved  and  original 
ZEM/ZEV  guidance  algorithms  and  to  compare  them  in 
every  aspect.  Simulation  results  show  that  although  the 
improved  ZEM/ZEV  guidance  law  has  slightly  less  landing 
precision  compared  to  the  original  one,  it  is  definitely  better 
in  avoiding  collision  to  the  planet  surface. 

Sliding  mode  control  (SMC)  is  a  hot  research  topic  in 
recent  years.  It  can  ensure  stability,  reliability  and  robust¬ 
ness  in  different  control  systems  under  actuator  faults, 
external  disturbances,  model  uncertainties  and  time-delay 
(Zhao  et  al.,  2014;  Wu  et  al.,  2012;  Liu  et  al.,  2011;  Lin 
et  al.,  2011;  Wu  et  al.,  2010;  Zhang  et  al.,  2010).  So  in 
Ebrahimi  et  al.  (2008),  Furfaro  et  al.  (2011,  2013),  the 
ZEM/ZEV  optimal  sliding-mode  guidance  law(OSG) 
which  is  a  combination  of  the  ZEM/ZEV  optimal  guidance 
law(OGL)  and  sliding-mode  theory  is  discussed.  However, 
from  massive  simulation,  the  OSG  does  not  show  any  sig¬ 
nificant  advantage  in  landing  precision  and  fuel- optimality 
compared  to  OGL  in  this  Mars  power  descent  problem. 
Besides,  the  ZEM/ZEV  optimal  sliding-mode  guidance 
law(OSG)  is  unpractical  due  to  the  chattering  problem 
which  is  common  in  sliding-mode  control.  So  we  will  not 
discuss  OSG  in  this  article. 

This  paper  is  organized  as  follows.  The  powered  descent 
problem  is  formulated  in  Section  2.  The  improved  optimal 
feedback  guidance  law  with  ZEM/ZEV  is  presented  in  Sec¬ 
tion  3.  Section  4  gives  simulation  study  and  some  discus¬ 
sions.  Finally,  some  concluding  remarks  are  given  in 
Section  5. 

2.  Problem  formulation 

2.7.  Equations  of  motion 

During  the  powered  descent  phase,  the  vehicle’s  altitude 
is  only  within  a  few  kilometers  and  its  velocity  is  also  very 
low.  So  the  aerodynamic  drag  and  lift  forces  can  be 
neglected.  The  equations  of  motion  during  the  powered 
descent  phase  for  a  Mars  lander  can  be  described  by 


r  =  v 

(1) 

v  =  g  +  a 

T 

(2) 

a  =  — 
m 

(3) 

T 

m  = - = 

nkTm  cos  6 

(4) 

c 

fpSe 

where  r  and  v  are  the  position  and  velocity  vectors,  a  is  the 
control  acceleration  provided  by  the  thrusters,  T  is  the 
thrusting  force  vector,  m  is  the  vehicle  mass,  and  g  is  the 
gravitational  acceleration  vector  acting  on  the  vehicle. 
During  the  powered  descent  phase,  g  is  considered  to  be 
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constant  due  to  the  little  variation  of  altitude.  These  vec¬ 
tors  are  all  3  x  1  column  vectors  expressed  in  a  nonrotating 
inertial  reference  frame  with  its  origin  at  the  target  landing 
site  of  Mars  (Guo  et  ah,  2013).  The  surface  fixed  reference 
frame  is  shown  in  Fig.  1.  It  is  shown  that  the  thrusting 
force  vector  T  can  be  expressed  as 


T  cos  a  sin  <7 

nkTm 

cos  0  cos  a  sin  o- 

T  =[TxTyTf  = 

T  cos  a  cos  g 

= 

nkTm 

cos# cos  a cos  fl¬ 

T  sina 

nkTm 

ees  Osin  a 

where  |T|  =  nkTmcos6  represents  the  two-norm  of  vector 
T,  which  is  the  magnitude  of  the  thrusting  force,  n  is  the 
number  of  identical  thrusters,  k  is  the  throttle  level  of  each 
thruster  and  Tm  is  the  maximum  available  thrust  magni¬ 
tude.  6  is  the  cant  angle  of  the  thrusters  to  net  thrust  direc¬ 
tion  (shown  in  Fig.  1).  a  is  the  thrust  angle  defined  as  the 
angle  between  the  thrusting  force  T  and  the  horizontal 
plane.  And  o  is  the  thrust  roll  angle  which  is  the  angle 
between  the  projection  of  T  on  the  horizontal  plane  and 
the  downrange  direction,  a  and  o  are  also  shown  in 
Fig.  1.  m  is  the  negative  mass  flow  rate  and  c  =  I spge  is 
the  engine  exhaust  velocity  with  Isp  and  ge  representing 
the  impulse  of  the  thrusters  and  earth  gravitational  acceler¬ 
ation  respectively. 

2.2.  Constraints 

During  the  powered  descent  phase,  the  lander  has  a  path 
constraints  with  the  boundary  conditions  as  follows  (Guo 
et  al.,  2013): 

r(?0)  =  r0  =  [x0  y0  zo]T)  r(tf )  =rf  =  [xf  yf  zf ]T  (6) 

vOo)  =  v0  =  [x0  y0  zo]T,  y(tf)  =\f  =  [xf  yf  z/]T  (7) 
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In  addition,  the  landing  guidance  law  should  also  ensure 
that  the  lander  will  stay  above  the  Martian  surface  all  the 
time  during  the  powered  descent  phase  in  order  to  avoid 
collision  with  the  surface  of  Mars,  which  is  a  precondition 
for  a  safe  and  successful  landing.  The  altitude  constraint 
can  be  expressed  as 

h  =  Ar  ^  0.  for  t0  ^  t  <  tf  (8) 

where  A  =  [0  0  1]. 

There  are  two  types  of  spacecraft  engines  -  power-lim¬ 
ited  and  thrust-limited  engines.  But  for  Mars  landing  mis¬ 
sion,  thrust-limited  engine  is  most  likely  to  be  used.  So  only 
thrust-limited  engine  will  be  considered  in  this  paper.  And 
the  engine  exhaust  velocity  of  the  thrust-limited  engine  is 
fixed  at  a  constant  value.  Theoretically,  the  throttle  level 
of  the  thrust-limited  engine  is  between  0  to  1,  but  in  real 
practice,  there  is  constraints  on  throttle  level.  In  this  paper, 
the  constraint  on  throttle  level  is  set  between  0.3  to  0.8, 
which  is  the  same  as  that  in  Acikmese  and  Ploen  (2007). 
So  the  thrust  magnitude  is  constrained  as 

0.3 Tm  ^  |T|  ^  0.8Tm  (9) 

where  Tm  is  the  maximum  available  thrust  magnitude.  So 
the  thrust  magnitude  constraint  will  form  a  control  satura¬ 
tion  problem  for  a. 

2.3.  Zero-effort  errors 

The  zero-effort  miss  at  time  t ,  ZEM(t),  is  the  distance 
(vector)  that  the  vehicle  will  miss  its  target  at  the  final  time 
tf  with  no  acceleration  command  exerting  on  the  vehicle 
after  the  current  time  t.  Similarly,  the  zero-effort  velocity 
error  at  time  t,  ZEV(t),  is  the  error  of  velocity  (vector)  at 
the  final  time  tf  if  no  acceleration  command  exerts  on  the 
vehicle  after  the  current  time  t.  Therefore, 
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ZEM(t)  =  Tf  -  r(/y),  a(r)  =  0,  for  t^T^tf  (10) 
ZEV(f)  =  \f  —  y{tf),  a(r)  =  0,  for  t  ^  t  ^  tf  (11) 

where  iy  and  V/  are  the  desired  final  position  and  velocity 
of  the  lander.  Considering  that  the  gravitational  accelera¬ 
tion  g  is  the  only  acceleration  exerting  on  the  vehicle  which 
is  also  assumed  to  be  a  constant  value,  the  ZEM(t)  and 
ZEV(t)  can  be  simplified  as  Guo  et  al.  (2013) 

ZEM  (0  =rf-(r  +  tgov  +  l-  t2g0  g)  (12) 

ZEV(t)  «  v7  -  (v  +  tgog)  (13) 

where  tgo  =  tf  —  t  is  the  time-to-go  before  reaching  the  ter¬ 
minal  state. 

3.  Improved  optimal  feedback  guidance  with  ZEM/ZEV 

For  optimization  problem,  the  classical  performance 
index  employed  in  Ebrahimi  et  al.  (2008),  Furfaro  et  al. 
(2011,  2013),  Guo  et  al.  (2013,  2011),  Hawkins  et  al. 
(2011,  2012),  Guo  et  al.  (2013)  is  as  follow 

1  ftf 

J  =  -  /  aTa  d t  (14) 

2  Jto 

The  Hamiltonian  function  is  defined  as 

H  =^a'a  +  pjv  +  pj(g  +  a)  (15) 

where  pr  and  py  are  the  costate  vectors  associated  with  the 
position  and  velocity  vectors,  respectively. 

In  Steinfeldt  et  al.  (2010)  and  D’Souza  and  An  (1997), 
an  additional  term  corresponding  to  the  final  time  tf  is 
added  to  the  classical  performance  index  (14).  This  repre¬ 
sents  a  trade-off  between  the  minimum  time  problem  and 
the  minimum  control  effort  problem. 

However,  the  classical  performance  index  and  Hamilto¬ 
nian  function  above  do  not  take  into  account  the  altitude 
constraints  in  (8)  which  may  cause  the  vehicle  to  crash 
the  ground.  In  order  to  rule  out  this  possibility,  we  propose 
a  novel  performance  index  which  add  a  term  corresponding 
to  the  altitude  constraints  in  (8)  to  the  classical  perfor¬ 
mance  index.  The  novel  performance  index  is  as  follow 


Remark  1.  Note  that  for  the  optimization  problem,  we 
want  to  minimize  the  performance  index,  also  called  the 
cost  function,  (14).  The  algorithm  we  will  use  to  minimize 
the  performance  index  with  the  help  of  Hamiltonian 
function  will  be  discussed  next.  So  if  we  add  a  term  —  Lir 
to  the  performance  index,  when  the  altitude  z  —  A r  is  below 
Zf,  the  performance  index  will  become  smaller  because  of 
an  extra  negative  term  —Li r.  Due  to  the  purpose  of  the 
algorithm  -  minimize  the  performance  index,  there  will  be 
an  extra  thrusting  force  acting  on  the  vehicle  automatically, 
which  will  make  it  to  fly  upwards,  to  make  the  extra  term 
—Li r  become  smaller,  until  the  altitude  reaches  Zf. 

As  the  gravitational  acceleration  g  is  assumed  to  be  the 
constant  value  in  this  problem,  the  optimal  control  equa¬ 
tion  and  costate  equations  are  as  follows 


dH 

=  o  a  =  -pc 

da. 

dH  1  T 

Pr  = 

-^  =  2^ 

dH 

Pr  = 

(19) 

(20) 
(21) 


The  solutions  of  the  optimal  control  equation  and  cost¬ 
ate  equations  are 

Pr=P r{tf)  (22) 

Pt,  =  -  \*ATiLg0  +  tgopr{tf)  +  p  ,(#/)  (23) 

3  =  -P„  =  \kAJ(o  ~  tgoVrhf)  -  P vhf)  (24) 

Integrating  the  acceleration  a,  the  velocity  and  position 
vectors  can  be  obtained  as 

V  =  ~Y2UTtlo  +  ploVrhf)  +  tgoVvhf)  ~  tgo g  +  V/  (25) 

r  =  4^8  ^T4>  riff)  v(tf)  +  \t2gog  -  tg0\f 

+ r/  (26) 

Combining  (25)  and  (26)  leads  to 


1  ftf  * 

J  =  -  /  (aTa  -  Li r) 

2  Jto 


At 


i  rtf . 

Jt0 

where  2  is  a  parameter  which  we  design  as  follow 

'0  z  >  Zf 


2  = 


c  z  <  z 


(16) 


(17) 


•f 


where  c  is  a  positive  constant.  The  Hamiltonian  function  is 
then  defined  as 


H  =  IaTa-tl4r  +  pJv  +  P:(g  +  a) 


(18) 


P  r(tf)  =  t2 


6(v/+v)  12(r  —  Yf)  1  T 


t3  +4AA  t» 


go  go 

2(2v/  +  v)  6(r  —  Tf)  1  X2 

24  > 


p vhf)  =  - v  2  —  +  g 

lgo  Lgo 


(27) 


(28) 


t  t2 

lgO  Lgi 

So  the  improved  optimal  guidance  law  can  be  obtained 


as 


_  6[rf  -  (r  +  ^v)]  2(\f  -  v)  1  T  2 

7L  1Z  g  +  24^  ^ 


(29) 


For  a  soft  landing  problem  with  vy  =  0,  the  improved 
OGL  has  the  following  form 


Please  cite  this  article  in  press  as:  Zhou,  L.,  Xia,  Y.  Improved  ZEM/ZEV  feedback  guidance  for  Mars  powered  descent  phase.  J.  Adv.  Space  Res. 
(2014),  http://dx.doi.Org/10.1016/j.asr.2014.08.011 


ARTICLE  IN  PRESS 


L.  Zhou,  Y.  Xia  /  Advances  in  Space  Research  xxx  (2014)  xxx-xxx 


5 


Table  1 

Initial  perturbations. 


Initial  perturbations 

Mean  value 

Standard  dev. 

Crossrange 

0  m 

50  m 

Downrange 

-2000  m 

150  m 

Altitude 

1500  m 

100  m 

Velocity  x 

0  m/s 

10  m/s 

Velocity  y 

100  m/s 

10  m/s 

Velocity  z 

-75  m/s 

5  m/s 

Mass 

1905  kg 

30  kg 

a_6(r/-r)  4v  l_  T  2 
tl„  t„  g  +  24^  go 


(30) 


Substituting  (12)  and  (13)  into  the  above  equation,  we 
get  the  improved  optimal  feedback  guidance  law  with 
ZEM/ZEV 


6  2 

a  =  ZEM - ZEV 


(31) 


Remark  2.  We  can  see  from  (29)— (31)  that  this  improved 
optimal  guidance  law  adds  an  extra  term  ^LiTt2go  to  the 
original  optimal  guidance  law  (OGL).  This  is  consistent 
with  our  design  purpose  of  performance  index  (16), 
because  when  the  altitude  z  ^  z/,  the  term  ^LiTt20 
becomes  nonzero  and  is  acted  as  an  extra  thrusting  force 
on  vehicle  which  helps  it  to  fly  upwards,  avoiding  collision 
to  the  planet  surface. 


Fig.  2.  Fuel  usage  for  various  tf. 

blow  the  dust  up  which  may  have  a  negative  impact  on  the 
vehicle.  The  navigation  system  on  the  lander  may  be  seri¬ 
ously  damaged  if  the  dust  is  blew  up  and  covers  the  cam¬ 
eras.  So  in  order  to  avoid  this  adverse  situation,  the 
powered  descent  phase  should  terminate  at  an  altitude 
above  the  Martian  surface.  The  Mars  Science  Laboratory 
(MSL)  mission  has  designed  a  sky  crane  landing  for  the 
final  touchdown  phase,  which  has  proven  to  be  a  good 
solution  to  this  problem. 

The  initial  perturbations  are  shown  in  Table  1  for 
Monte  Carlo  simulations. 


4.  Simulation  results 

4.1.  Parameter  design 

In  this  paper,  the  parameters  of  the  vehicle  for  simula¬ 
tion  are  mainly  from  Acikmese  and  Ploen  (2007),  which 


are  as  follows 

Wo  =  1905  kg,  g  =  [0  0  -  3.71 14]T  m/s2 
n  =  6.  6  =  27  deg  (32) 

Isp  =  225  s,  ge  =  9.807,  Tm  =  3100  N 

The  initial  position  and  velocity  of  the  lander  are  given 
as  follows  (Acikmese  and  Ploen,  2007) 

r(t0)  =  r0  =  [0  —  2000  1500]T  m  (33) 

v(t0)  =  Vo  =  [0  100  -  75]t  m/s  (34) 

In  this  paper,  the  target  terminal  position  and  velocity  are 
set  as 

r {tf)  =  iy  =  [0  0  100]T  m  (35) 

y{tf)  =  \f  =  [0  0  0]T  m/s  (36) 


The  target  terminal  altitude  is  set  to  100  m  for  the  final 
touchdown  phase  to  execute.  If  the  target  terminal  position 
is  set  to  the  surface  of  Mars,  the  air  flow  of  the  engines  will 


4.2.  Determination  of  optimal  time-to-go 

From  (31),  it  is  shown  that  the  time-to-go  tgo  =  tf  —  t 
has  an  important  effect  on  the  performance  of  the  optimal 
feedback  guidance  laws  (original  and  improved  OGL). 
Unless  the  total  time-to-go  of  the  powered  descent  phase 


Fig.  3.  Minimum  altitude  for  various  tf. 
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Table  2 

Landing  precision  for  various  tf. 

L.  Zhou, 

Y  Xia  /  Advances  in  Space  Research  xxx 

(2014)  xxx-xxx 

Mean  value  vs  tf 

30  s 

35  s 

40  s 

45  s 

Crossrange 

0.1292  m 

-0.0247  m 

-1.0710  x  10-4  m 

4.1191  x  10-6  m 

Downrange 

58.9524  m 

1.7807  m 

0.0149  m 

1.1241  x  10-4  m 

Altitude 

23.3845  m 

97.5814  m 

99.9726  m 

99.9997  m 

Velocity  x 

0.0311  m/s 

-0.0105  m/s 

-6.3293  x  10~4  m 

3.2170  x  10~4  m/s 

Velocity  y 

13.9600  m/s 

1.0585  m/s 

0.0346  m/s 

—0.0074  m/s 

Velocity  z 

-16.5214  m/s 

-1.4892  m/s 

-0.0728  m/s 

—0.0042  m/s 

Mass 

1742.4  kg 

1718.7  kg 

1703.9  kg 

1690.8  kg 

Table  3 

Landing  statistics  with  OGL. 


Landing  statistics 

Mean  value 

Standard  dev. 

Crossrange 

5.9124  x  10-6  m 

8.0346  x  10-5  m 

Downrange 

1.3202  x  10“4m 

1.4106  x  10“4  m 

Altitude 

99.9998  m 

3.9655  x  1CT4  m 

Velocity  x 

5.8410  x  10“4  m/s 

0.0049m/s 

Velocity  y 

-0.0071  m/s 

0.0148  m/s 

Velocity  z 

—0.0014  m/s 

0.0330  m/s 

Mass 

1693.1  kg 

28.4988  kg 

Table  4 

Landing  statistics  with  improved  OGL. 

Landing  statistics 

Mean  value 

Standard  dev. 

Crossrange 

1.7195  x  10“5m 

2.3863  x  10"4  m 

Downrange 

1.0136  x  10“4m 

5.3703  x  10"4  m 

Altitude 

99.9995  m 

0.0027  m 

Velocity  x 

4.2053  x  10"4  m/s 

0.0057  m/s 

Velocity  y 

—0.0070  m/s 

0.0120  m/s 

Velocity  z 

—0.0041  m/s 

0.0618  m/s 

Mass 

1691.3  kg 

27.0150  kg 

is  constrained,  the  time-to-go  tgo  in  every  guidance  cycle 
can  be  determined  to  achieve  the  optimal  performance  in 
fuel  consumption  and  precision  landing.  But  this  is  not 
easy  considering  the  large  amount  of  time  which  will  be 
spent  on  finding  the  optimal  time-to-go.  Guo  et  al.  (2011) 
proposed  a  numerical  method  to  calculate  the  optimal 
time-to-go  for  the  power-limited  engine.  But  for  a  thrust- 
limited  engine,  minimizing  J  (see  (14))  does  not  mean 
minimizing  the  fuel  usage,  so  the  result  obtained  by  the 
numerical  method  mentioned  above  is  not  the  optimal 
time-to-go  for  the  thrust-limited  engine.  As  a  result,  the 
optimal  ZEM/ZEV  guidance  laws  discussed  in  this  paper 
are  not  fuel-optimal  in  nature,  but  the  fuel  usage  can  be 
regarded  as  the  performance  measure  of  interest  (Guo 
et  al.,  2011).  In  this  section,  the  impact  of  time-to-go  on 
the  performance  of  the  optimal  guidance  laws  will  be  dis¬ 
cussed  and  the  optimal  time-to-go  of  the  Mars  powered 
descent  phase  will  be  determined  by  trial.  Note  that  in  the 
following  we  will  use  the  simulation  results  of  the  original 
OGL  to  study  and  determine  the  optimal  time-to-go. 

In  order  to  find  the  optimal  time-to-go  which  can  mini¬ 
mize  the  fuel  usage,  the  fuel  usage  for  various  tf  is  shown  in 
Fig.  2  below.  From  Fig.  2,  it  is  shown  that  the  fuel  usage  is 


500 


-2500  _ i _ i _ i _ ■  _ i _ i _ i _ 

-250  -200  -150  -100  -50  0  50  100  150  200 


Crossrange  <m) 

Fig.  4.  Trajectories  of  OGL. 

positive  proportional  to  tf.  So  the  less  time-to-go  is,  the  less 
fuel  will  be  consumed.  But  the  optimal  time-to-go  cannot 
be  too  small  considering  the  constraints  on  control  forces. 
Besides,  the  no-subsurface  constraint  will  also  be  used  to 
find  the  optimal  time-to-go.  The  minimum  altitude  in  every 
trajectory  for  various  tf  is  shown  in  Fig.  3.  From  Fig.  3,  we 
can  see  that  when  the  final  time  tf  is  greater  than  57  s,  the 
minimum  altitude  will  be  below  100  m,  which  is  undesir- 

500  - 1 - T - 1 - 1 - 


-2500  _ i _ i _ i _ i _ i _ i _ i _ i _ 

-200  -150  -100  -50  0  50  100  150  200  250 

Crossrange  (m) 

Fig.  5.  Trajectories  of  improved  OGL. 
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Fig.  6.  Landing  dispersion  of  OGL. 
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Fig.  9.  Altitude  vs  time  of  improved  OGL. 
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Fig.  7.  Landing  dispersion  of  improved  OGL. 


Fig.  10.  Vertical  velocity  vs  time  of  OGL. 


able  considering  the  greater  chance  of  collision  if  the  target 
site  is  at  the  Martian  surface.  So  the  final  time  tf  is  better  to 
set  below  57  s. 

Another  problem  that  needs  to  be  concerned  is  the  pre¬ 
cision  landing  requirement.  If  the  time-to-go  is  set  too  low, 
the  thrusting  force  will  saturate,  which  may  not  satisfy  the 
precision  landing  requirement.  The  landing  precision  for 
various  tf  with  initial  perturbations(see  Table  1)  will  be 
shown  in  Table  2  below. 

The  data  above  shows  that  when  the  final  time  tf  is 
below  40  s,  the  final  velocity  is  greater  than  1  m/s,  which 
cannot  satisfy  the  soft  landing  requirement.  When  the  final 
time  tf  is  greater  than  40  s,  the  landing  precision  is  good 
enough  to  satisfy  the  requirement.  The  position  and  veloc¬ 
ity  errors  have  a  magnitude  at  least  of  10-3.  Besides,  the 
final  mass  for  various  tf  agrees  with  the  tendency  shown 
in  Fig.  2  that  the  fuel  usage  is  direct  proportional  to  the 
final  time  tf.  So  in  order  to  be  fuel-optimal,  the  final  time 
should  be  set  between  40  s  to  45  s. 


Fig.  8.  Altitude  vs  time  of  OGL. 
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Fig.  1 1 .  Vertical  velocity  vs  time  of  improved  OGL. 


Fig.  14.  Thrust  angle  vs  time  of  OGL. 
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Fig.  12.  Commanded  acceleration  vs  time  of  OGL. 
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Fig.  15.  Thrust  angle  vs  time  of  improved  OGL. 
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Fig.  13.  Commanded  acceleration  vs  time  of  improved  OGL. 
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Fig.  16.  Thrust  roll  angle  vs  time  of  OGL. 
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The  main  difference  of  the  original  and  improved  OGL 
lies  in  the  altitude  profiles.  Fig.  8  shows  that  about  one  fifth 
of  the  cases  have  a  minimum  altitude  below  100  m  during 
the  power  descent  phase.  In  some  cases,  the  minimum  alti¬ 
tude  reaches  0  m,  or  even  below  0  m.  If  we  set  the  final  alti¬ 
tude  below  100  m,  let’s  say  0  m  -  planet  surface,  in  many 
cases  the  vehicle  will  crash  the  ground  before  reaching 
the  target  landing  place.  However,  we  can  see  from  Fig.  9 
that  only  several  cases  have  a  minimum  altitude  below 
100  m  during  power  descent  and  they  are  just  slightly 
below  100  m.  In  most  of  the  cases,  the  minimum  altitude 
is  just  100  m,  which  forms  a  horizontal  straight  line  on 
100  m  as  shown  in  Fig.  9.  So  this  demonstrates  that  the 
improved  optimal  guidance  law  is  effective  in  preventing 
subsurface  flight,  in  other  words,  collision  to  the  surface. 


5.  Conclusions 


0  5  10  15  20  25  20  35  40  45 

time  (s) 

Fig.  17.  Thrust  roll  angle  vs  time  of  improved  OGL. 
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The  final  time  tf  will  be  set  to  45  s  for  the  Monte  Carlo 
simulation  in  the  next  subsection  to  leave  some  control 
margin  for  the  assumed  initial  perturbations. 

4.3.  Monte  Carlo  simulations 

In  order  to  assess  the  robustness  and  reliability  of  the 
original  and  improved  OGL  and  to  compare  them  in  every 
aspect,  Monte  Carlos  with  300-dispersed  cases,  varying  ini¬ 
tial  states  and  model  errors,  were  run  using  Matlab  simula¬ 
tion  for  the  two  guidance  algorithms.  The  guidance  cycle  is 
set  to  10  ms.  The  initial  dispersions  are  listed  in  Table  1  for 
both  two  algorithms. 

In  addition,  errors  in  mass,  thrust,  gravitational  acceler¬ 
ation  and  wind  influence  will  also  be  added  in  simulation  as 
a  perturbation  p  on  the  acceleration,  p  is  assumed  to  be 
p(t)  =  Psin(f  f),  where  P  =  0.2a. 

The  landing  statistics  of  the  original  OGL  and  improved 
OGL  are  shown  in  Tables  3  and  4.  And  the  comparison  of 
the  original  and  improved  OGL  on  trajectories,  landing 
dispersions,  altitude,  vertical  velocity,  commanded  acceler¬ 
ation,  thrust  angle  and  thrust  roll  angle  are  shown  in 
Figs.  4-17. 

As  is  shown  in  Tables  3  and  4,  although  the  landing  pre¬ 
cision  of  the  original  OGL  is  slightly  better  than  that  of  the 
improved  one,  it  is  neglectable  considering  the  large  size  of 
the  lander.  So  we  can  say  that  the  original  and  improved 
OGL  have  almost  identical  performance  in  precision  land¬ 
ing  and  fuel  consumption.  And  the  profiles  of  altitude,  ver¬ 
tical  velocity,  commanded  acceleration,  thrust  angle  and 
thrust  roll  angle  are  nearly  the  same.  This  is  reasonable 
considering  the  nature  of  the  improved  OGL  proposed  in 
this  paper.  Only  when  the  altitude  is  below  the  desired  final 
altitude,  an  additional  term  is  added  to  the  original  optimal 
guidance  law  to  let  it  fly  upwards,  avoiding  collision  to  the 
ground.  When  the  altitude  is  above  the  desired  final  alti¬ 
tude,  the  improved  OGL  is  identical  to  the  original  one. 


For  Mars  power  descent  problem,  an  improved  optimal 
guidance  law  with  ZEM/ZEV  is  proposed  in  this  paper  to 
address  the  altitude  constraints.  This  method  can  be  fuel- 
optimal,  at  the  same  time  satisfying  the  precision  landing 
and  no-subsurface  flight  requirement.  Besides,  the  optimal 
time-to-go  has  been  studied  and  determined  by  simulation 
method.  Simulation  results  show  that  both  the  original  and 
improved  optimal  guidance  law  have  robust  and  reliable 
performances  under  perturbations.  And  simulation  also 
demonstrates  the  effectiveness  of  the  improved  OGL  in 
avoiding  collision  to  the  planet  surface. 
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